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Abstract 

We describe two algorithms for finding representatives of the nilpotent orbits of 
a 0-group. The algorithms have been implemented in the computer algebra sys- 
tem GAP (inside the package SLA). We comment on their performance. We apply 
the algorithms to study the nilpotent orbits of ^-groups, where 9 is an N-regular 
automorphism of a simple Lie algebra of exceptional type. 

1 Introduction 
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^ ■ Let g be a simple complex Lie algebra, and let G be a connected algebraic group with Lie 
0> ! algebra q. Then G acts on g. The orbit structure of the action of G on g has been studied 
^ ' in detail (we refer to [6] for an overview). In particular, the nilpotent orbits have been 
classified, using a correspondence between nilpotent G-orbits, and G-conjugacy classes of 
^ ' s[2-triples. 

■ Kostant and Rallis (tl9j) considered decompositions of the form g = go©gi, where g^ is 

the eigenspace of an involution 9 of g, corresponding to the eigenvalue (—1)*. Let Go C G 
be the connected subgroup with Lie algebra go. Then the reductive group Go acts on gi, 
and again the question arises as to what its orbits are. In [19], among many other things, 
it is shown that there are finitely many nilpotent orbits, and a correspondence between 
nilpotent orbits and s[2-triples, analogous to the one for g, is established. 

This was generalised by Vinberg in the 70 's ([2S], [SD])- He considered the decomposition 
of g relative to an automorphism 6 of order m (or relative to a 1-parameter group of 
automorphisms if m = cxd). Here 

= 0O©01©---©0m-l, 

where g^ is the eigenspace of 6 corresponding to the eigenvalue u^, where a; is a primitive 
m-th root of unity. Again we get a reductive connected algebraic group Go, with Lie 
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algebra go, acting on gi. The group Gq, together with its action on gi, is called a 6'-group. 
In the sequel we will, by a slight abuse of language, also call Gq a 6'-group; the action on 
01 is always understood. 

Example 1 Let g = sl^lC); then G = SL4(C) acts on g by conjugation. Let u; G C be 
a primitive third root of unity. Let 6 be the automorphism of order 3 of g, given by the 
following matrix 
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Here, if on position (i, j) there is uj^, then 0{eij) = uj'^Cij, where ejj is the matrix with 
a 1 on position and zeroes elsewhere. Let hi = Ci^i — Cj+i^j+i. Then we see that 

00 is spanned by hi, /12, h^, 63^4, 64^3; which means that 0o = sl2{C) © T2 (where T2 
denotes the subalgebra spanned by hi,h2). Furthermore, 0i is spanned by 61^2, ^2,3, 62,4, 
63,1 and 64,1. As go-module (and hence as Go-module) gi splits as a direct sum of two 2- 
dimensional modules (spanned respectively by 62,3,62,4 and 63,1,64,1) and a 1-dimensional 
module (spanned by 61,2). 

It is of interest to study 6'-groups for a number of reasons. Firstly, they form a class 
of algebraic groups for which it is not completely hopeless to list the orbits. Secondly, 
many interesting representations of algebraic groups arise as 6'-groups (for example SL(9, C) 
acting on A^(C^), which is studied in [31j). 

An orbit G^e (with 6 G 0i) is called nilpotent if is contained in its closure. This 
happens if and only if e is nilpotent as an element of (that is, the adjoint map adg(6) 
is nilpotent). The results of Vinberg show that also here there are a finite number of 
nilpotent Go-orbits in gi. Secondly, there is a correspondence between nilpotent orbits and 
s[2-triples. Moreover, in [28], [30], Vinberg developed a method for obtaining the nilpotent 
Go-orbits in 0i. It is the objective of this paper to describe algorithmic methods, that can 
be implemented on computer, for this purpose. 

The conjugacy classes (in Aut(g)) of the finite order automorphisms of g have been 
classified by Kac ([IB], see also [TB]) in terms of so-called Kac diagrams. The connected 
component of the identity of Aut(0) is the group of inner automorphisms of g. It is 
generated by the exp(adg(6)), for 6 G g nilpotent. We have that an automorphism 9 is 
inner if and only if go contains a Cartan subalgebra of g (this follows, for example, directly 
from [16j, Chapter X, Theorem 5.15). 

As said before, the purpose of this paper is to describe algorithms for listing the nilpo- 
tent orbits of a ^-group. For this, our computational setup is as in [12]. In particular, we 
assume that g is given by a multiplication table relative to a Chevalley basis. Then all 
structure constants in the multiplication table are integers. From the Kac diagram of an 
automorphism 6 it is straightforward to compute the matrix of 6 relative to the given basis 
of g (this follows directly from the description of the finite order automorphisms given in 
[TB] , Chapter X, Theorem 5.15). Then by linear algebra we can construct bases of go and 
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01. In this paper we will focus on the case where 6 is an inner automorphism of g. This is 
by far the most common case. Secondly, in this case a Cartan subalgebra of go will also be 
a Cartan subalgebra of q. This implies that go and gi are generated by root spaces of g. 
So, although the matrix of 9 has coefficients in Q(ti;), the spaces go, gi are defined over Q. 
Therefore, all computations are completely rational, i.e., the field elements that appear all 
lie in Q. 

Several listings of nilpotent orbits of ^-groups have appeared in the literature. We 
mention pT] (g of type Eg, the order of 9 equal to 3, go of type Ag), p] (g of type Eg, the 
order of 9 equal to 2, go of type Ds), P (g of type Es, the order of 9 equal to 5, go of type 
2^4 and g of type E^, the order of 9 equal to 3, go of type A2 + A^), [26j (g of type E^, the 
order of 9 equal to 4, go of type Ai + 2A3). (We remark that in those references all orbits 
are classified, so also the semisimple orbits, and the ones of mixed type). Furthermore, in 
[22] and [23] the nilpotent orbits are listed for the cases where g is of exceptional type, and 
the order of ^ is 2. 

Littelmann ([20]) has devised an algorithm to obtain nilpotent orbits of 6'-groups. His 
algorithm is tailored towards the case where the 6'-group is defined by a Z-grading. We note 
that for Z-graded algebras g, a Cartan subalgebra of go will also be a Cartan subalgebra of 
g ([30], §1-4). For this reason, the algorithms in this paper work unchanged for Z-gradings; 
however, we will not explicitly consider these. The main step of Littelmann's algorithm 
consists of computing sets w{A), where A is a basis of the root system, and w runs through 
a set of representatives the right cosets of a Weyl subgroup Wq of the Weyl group W. For 
this reason, the algorithm will behave in a similar manner to our first method (described 
in Section [3]). That is, it will work well when the index of Wq in W is small. 

In [27], Popov has developed an algorithm for computing the Hesselink strata of the 
nuUcone of a representation of a reductive algebraic group. If the group in question is 
a ^^-group then this yields an algorithm for computing the nilpotent orbits. Experiments 
with an implementation of the algorithm (due to A'Campo and Popov) however suggest 
that for the special case of 6'-groups the more specialised methods of this paper perform 
better. (For an example, let 9 be the N-regular inner automorphism (see Section [6]) of 
order 2 of the Lie algebra of type Eq. Then go is of type A^ + Ai, and the module gi is 
irreducible with highest weight (0,0,1,0,0,1). For this case Popov's algorithm used 831 
seconds, wheras the methods of Sections [3] and H] needed respectively 5 and 93 seconds.) 

This paper is organised as follows. In Section [2] we describe two algorithms that we 
need for tasks related to Weyl subgroups of Weyl groups. In Section [3] we describe the first 
method for listing nilpotent orbits of ^-groups. It uses the classification of the nilpotent 
orbits of g. For each such orbit it is decided whether it intersects with gi, and if this is 
the case, the Go-orbits of the intersection are determined. This method works well if the 
index of the Weyl group of go in the one of g is small; and the method ceases to work well 
if this index is large. In that case, the dimensions of the spaces go and gi will be small 
(compared to the dimension of g). In Section H] we describe a second method, which is 
based on Vinberg's theory of carrier algebras. This method works well if the spaces go 
and Qi are small-dimensional. In this sense the two algorithms complement each other. 
All algorithms have been implemented in the computer algebra system GAP4 ([ID]), inside 
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the package SLA ([I3]). In Section [5] we briefly report on practical experience^ with these 
implementations. In Section [6l we apply the algorithms to study the nilpotent orbits in 
case 9 is an N-regular automorphism of g, where g is of exceptional type, of orders between 
2 and the Coxeter number. 

Acknowledgements: I would like to thank Alexander Elashvili and Oksana Yakimova 
for inspiring conversations on the subject of this paper. 

2 Algorithms for Weyl subgroups 

Let Let $ be a root system, with basis of simple roots A = {ai, ...,«;}. Then the group 
generated by the reflections Si = Sa^ is the Weyl group of $, which we denote by W. Let 

. . . , /3s be roots in $ that form a basis of a root subsystem of $. Then the Weyl group 
Wq of \E' is a subgroup of W, generated by the reflections Sj3^. The group Wq is said to be a 
Weyl subgroup of W. In this section we describe algorithms for two tasks concerning Weyl 
subgroups of Weyl groups: finding a set of representatives of the right cosets of Wq in W, 
and checking whether two sets of weights are conjugate under Wq. This is the content of 
the first two subsections. In the third subsection we describe an application to listing the 
TT-systems in up to ly-conjugacy. 

It is straightforward to see that the root system has a basis consisting of positive 
roots in $. Therefore, in the sequel we assume that the are positive roots. 

Let P denote the weight lattice of $, spanned by the fundamental weights Ai, . . . , A;. 
Let ( , ) be a VT-invariant inner product on the space spanned by P. For two weights A, fi 
we set (A,/i^) = 2(A, Let F denote the set of weights fi E P with (/x, a^) > 0. 
This set is called the fundamental Weyl chamber of W. It is known that every IV-orbit 
has a unique point in F. The fundamental Weyl chamber Fq of Wq is the set of all fi E P 
with > 0, for all i. Every Wo-orbit in P contains a unique point in Fq. 

For a w E W we denote by £{w) the length of w, that is, the length of a reduced 
expression for w as a word in the generators Sj. Similarly, for u G Wq, io{u) will denote 
the length of a reduced expression of u in the generators s^-. 

2.1 Listing right cosets 

It is known (see ^) that every right coset of a Weyl subgroup of a Weyl group has a 
unique representative of shortest length. Also there are algorithms known for finding a 
set of shortest coset representatives. (For example the computer algebra system Magma 
([3j), and the GAPS package Chevie ([H]) contain implementations of such an algorithm.) 
However, I have not been able to find a reference for such an algorithm in the literature. 
So for reasons of completeness, this section contains a characterisation of these shortest 
representatives, that also yields an algorithm to find them. 

^AU timings reported in this paper have been obtained on a 2GHz machine, with 1GB of memory for 
GAP. 
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Lemma 2 Let w & W be such that w~^{(3i) > for 1 < i < s. Then w is the unique 
element of shortest length in the coset WqW. Moreover, every coset WqU contains an 
element with this property. 

Proof. We claim the following: let v G W, and let P be one of the Pi. Suppose that 
v-\p) > 0. Then £{spv) > £{v). 

First we note that this is equivalent to ilv'^sp) > i{v~^). We use the fact that for 
u E W we have that i{u) is equal to the number of positive roots a E ^'^ such that 
u{a) < 0. Let 

[/ = {a e I sp{a) e $+}. 

Then permutes U, hence the number oi a E U with v''^{a) < is equal to the number 
of a in f/ with v~^Sf3{a) < 0. So and v~^sp receive the same "contribution" from U 
towards their lengths. Set 

1/ = {a e I S/3{a) G 

Then — permutes V. This implies that the number of a E V with v^^Sf3{a) < is equal 
to the number of a E V with v~^{a) > 0. We show that this last number is strictly bigger 
than the number of a eV with v~^{a) < 0. This then imphes the claim. 

Set Ml ^ {a eV \ {fi, a) < 0}, and M2 ^ {a E V \ (//, a) > 0}. Let a E Mi. Prom 
sp{a) < we get that (q;,/3^) > 0. But then 

v-\sp{a))=v-\a)-{a,P'')v-\p). 

Now since v~^{a) < and v~^{j3) > we get that v^^sp{a) < 0. It follows that —5^3 maps 
Ml into M2. But j3 G M2 does not lie in — S/3(Mi). Hence IM2I > |Mi| and we are done. 

Now suppose that w~^{Pi) > for 1 < -i < s. Let u E Wq be such that i{uw) > i{w). 
Again let P be one of the with £o{spu) > ioiu). The latter condition implies that 
u'^{P) is a linear combination of the Pi with non-negative integer coefficients. Hence 
w^-^u^-^i^P) > 0. Therefore, by the claim above, i{s0uw) > £{uw). Now by induction on 
io{u) wc sec that £{uw) > £{w) for all u E Wq, u ^ 1. 

Let V E W, and P one of the Pi. If v~^{P) < 0, then by arguments similar to the ones 
used before, we prove that i{spv) < i{v) (this time —sp maps M2 into Mi). Hence, if 
u E Wqw does not have the property that u{Pi) > for 1 < i < s, then we can find Pi with 
£{si3^u) < £{u). Continuing, we eventually find an element in the coset of minimal length. 

□ 

We say that a w G is a shortest representative if it is the unique representative of 
shortest length of the coset Wqw. 

Lemma 3 Let w be a shortest representative. Write w — w'sa^, where £{w') — £{w) — 1. 
Then w' is a shortest representative. 

Proof. If not, then we can write w' = Wiw", with Wi E Wq and £{w") < £{w'). Hence 
w and w"sa^ lie in the same right I^o-coset. Hence w"sa^ is not a shortest representative. 
Therefore we can write w"sai = W2w"', with W2 E Wq and w'" a shortest representative. 
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£{w"') < i{w"saj. But then w = wiW2w"', and this imphes w = w'". But i{w"') < i{w") < 
£{w') < i{w), which is a contradiction. □ 
Now let Rk denote the set of shortest representatives of length k. Lemmas [2] and [3] lead 
to the following algorithm for computing Rk+i from Rk- Initially we put -R^+i = 0- Then 
for 1 < z < / and w E Rk we do the following: if £(wSq,-) > i{w) and wSaXPj) > for 
1 < j < -5 then add wSq,. to Rk+i- 

Remark. The implementation of this algorithm in GAP4 appears to work well. For ex- 
ample for Wq of type 2A4 inside W of type Eg, (this is obtained by taking the set of simple 
roots, adding the lowest root, and deleting an appropriate simple root), the GAP4 imple- 
mentation takes 2.1 seconds, whereas Chevie and Magma V2. 14-11 need respectively 9.4 
and 89.7 seconds. In this example there are 48384 cosets. 



2.2 Checking conjugacy 

Let Fi = {/ii, . . . , fim} and F2 = {Ai, . . . , Am} be two subsets of P. In this section we 
describe how we can check efficiently whether there exists a w E Wq with ^^(Fi) = F2. For 
this we first focus on the problem of deciding whether there is a w E Wq with w{fii) = Aj 
for 1 < i < m. 

First we observe that it is straightforward to compute, for given ft E P, a w E Wq and 
a A G Fo with w{fi) = A. (Note that A is uniquely determined by fi.) Indeed, we find 
the smallest index i with {fi,Pi) < 0. If there is no such i then ft E Fq and we are done. 
Otherwise, we set /ii = S/3-(/i), and we continue with /ii in place of fi. This algorithm 
terminates as /ii > /i in the usual partial order on P (which is defined by z/ < 77 if 77 — i/ is 
a sum of positive roots). Furthermore, by tracing the sp. that we applied we find w. 

This means that the problem is easily solved if the sets Fj have only one element. 
Indeed, we compute u,v E Wq such that u(yUi) and v{Xi) lie in Fq. If u(yUi) = f(Ai) 
then they are conjugate, and w = v~^u is such that w(yUi) = Ai. Otherwise they are not 
conjugate. 

If the Fj are larger, then, as a first step, we decide whether there is a wi E Wq with 
wi(/ii) = Ai. If there is no such wi then our problem has no solution, and we stop. 
Otherwise we compute such a wi. 

Now we need an intermezzo on stabilisers. Let v E Wq be such that f (Ai) = X E Fq. 
For a weight u we consider its stabiliser Stabiyo('^) = {u E Wq \ uiy) = z/}. We note that 
we have an isomorphism a : Stabvi/o(Ai) — > Stabpi/Q(A), by a{u) = vuv'^. Now as A G Fq 
we have that Stabvi/o(A) is generated by the s^^ with (A,/5/) = (cf. |12J, Proposition 
8.3.9 - there it is proved for the full Weyl group, but the proof goes through also for Weyl 
subgroups). Hence StabvKo(Ai) is generated by all v~^Si3^v = where is such that 

(A,/5/) = 0. We conclude that Stabvi/o(Ai) is a Weyl subgroup of W; moreover, we can 
compute the refiections that generate it. 

It is straightforward to see that Stabiyo(Ai)wi is exactly the set of elements of Wq that 
send fii to Ai. Now set /i^ = wi(/ij) for i > 1. Then by induction on the size of Fj we can 
decide whether there exists a w E Stabvi/o(Ai) with w{fi'^ = Aj for 2 < i < m. Now if such 
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a w exists, then v = wwi G Wq has the property that f (/ij) = Aj. Otherwise, such a v does 
not exist. 

Now we return to the more general problem, i.e., to decide whether there exists a 
w G Wo with w{Ti) = T2- We assume that the fii and Aj are ordered in such a way that 
the matrix Bi = {{fii, /Xj))^-^^ is equal to the matrix B2 = ((Aj, Aj))™-^]^. Then we compute 
all permutations r of Fi that leave Bi invariant, i.e., such that 

= {T{fii),T{fij)) for l<i,j <m. 

Then for each such r we check whether there is w G Wq with wT{iii) = Aj. 
Remark. This algorithm works rather well in practice. First of all, usually there are 
not many permutations r that leave Bi invariant. Secondly, the basic operation of the 
algorithm is to compute a A G Fq conjugate to a given weight; and this can be done in 
few steps. We have used an implementation in GAP4 of Dynkin's algorithm for classifying 
so-called vr-systems (cf. Section [2l3i) . up to VT-conjugacy, in the Lie algebra of type E^. 
A list of 76 vr-systems was constructed, and the algorithm for deciding conjugacy under 
the Weyl group of type Es was called 3873 times. The total time used was about 71 seconds. 

2.3 Dynkin's 7r-systems 

Let F C $; then F is called a vr-system if 
CI) for all a,P eT we have a - /? ^ $, 
C2) F is linearly independent. 

We have that F C $ is a basis of a root subsystem of $ if and only if it is a vr-system. 

In [8] Dynkin gave a neat algorithm to classify vr-systems of maximal rank (i.e., of rank 
equal to the rank of $). This works as follows. Let F be a vr-system, and D G T a subset 
corresponding to a connected component of the Dynkin diagram of F. Then D is a basis 
of a root subsystem of $. To D we add the lowest root of that root subsystem. Secondly, 
we erase a root from D, different from the one added. This yields a vr-system F'; which 
is said to be obtained from F by an elementary transformation. Dynkin showed that all 
vr-systems of maximal rank can be obtained (up to VT-conjugacy) from A, by a series of 
elementary transformations. So in order to get a list of all vr-systems, up to W^-conjugacy, 
one does the following: 

1. Find all maximal vr-systems that can be obtained from A by perfoming elementary 
transformations . 

2. From this set erase VT-conjugate copies, to obtain the set M. 

3. Let M' be the set obtained from M by adding all subsets of each element of M. 

4. From M' erase VT-conjugate copies. 
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We remark that for checking whether two vr-systems are H^-conjugate, the algorithm 
from Section 12.21 can be used. 

Remark. Using our implementation in GAP of this algorithm, we have obtained the same 
tables for the root systems of exceptional type as Dynkin (|8j). For the root systems of 
classical type this algorithm has been applied by Lorente and Gruber (|21])- However, 
with our implementation of the algorithm we obtained tables that for many root systems 
contain more 7r-systems. 



3 Listing nilpotent orbits I 

Now we return to the set up of Section [TJ That is, 6 is an automorphism of q of order m, 
and g = ©^o^fl* the corresponding Z/mZ-grading. We describe a method for finding 
representatives of the nilpotent Go-orbits in gi. 

3.1 Preliminary lemmas 

We start with three lemmas which, for the case m = 2, have been proved in [19]. The 
proofs for general m are entirely similar, and therefore we omit them. A triple e, /) of 
elements of g, with [h, e] = 2e, [h, /] = —2/, [e, /] = h, is called an s[2-triple. 

Lemma 4 Let e G gi he nilpotent. Then there are h E Qq and f G g_i such that {h, e, /) 
is an 3i2-triple. 

The proof is the same as the first part of the proof of [19], Proposition 4. This lemma is 
also part of the content of [30], Theorem 1. We call an s[2-triple {h, e, /) with the properties 
of LemmaHla norma/ s[2-triple, in analogy to [19J. The group Go acts on normal s[2-triples 
by g ■ {h,ej) = {g ■ h,g ■ e,g ■ /). 

Lemma 5 Let A be the set of nilpotent Go-orbits in Qi. Let B be the set of Gq- orbits of 
normal si2-triples. Let ip : A ^ B be defined as follows: ^{Gq ■ e) = Gq ■ {h,e,f), where 
{h,e,f) is any normal 5l2-triple containing e. Then (f is well-defined, and bijective. 

The proof is analogous to the second part of the proof of [19j, Proposition 4. 

Lemma 6 Let {h,e,f) and (/ii,ei,/i) be normal 5l2-triples. They are conjugate under Gq 
if and only if h and hi are. 

Here the proof follows the one of [19], Lemma 4. 

We say that an /i G go is normal if it lies in a normal 512-triple {h, e, /). By Lemmas HI 
[5] and El listing the nilpotent Go-orbits in gi is equivalent to listing the Go-orbits of normal 
h G go. 
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3.2 Deciding normality 



We call a.n h E Qo admissible if there are e, f E Q such that {h, e, /) is an s[2-triple. Let f)o 
be a fixed Cartan subalgebra of go- In this section we describe an algorithm for deciding 
whether a given admissible /i G go is normal. 

Proposition 7 Let h E Qq be admissible and write Qi{k) = {x G 0j | [h,x] = kx}. Set 

U = {ueQ,{2) I [go(O),n]=0i(2)}. 

Then U is dense in 0i(2). Moreover, an e (zU lies in a normal sl2-triple {h, e, /) or there 
is no normal sl2-triple containing h. 

Proof. Consider the subalgebra 

^ ~ Qi mod m(2i)- 

Then a is reductive. Indeed, let k be the Killing form of g. By standard arguments (i.e., 
analogous to [17], Chapter IV, §1.1) we see that t^idiik), gj{l)) = unless j = —i and 
/ = —k. So, since g is the direct sum of the gi{k), the restriction of k to gi{2i) © 2i) 
is nondegenerate. This implies that the restriction of k to a is nondegenerate, and hence 
that a is reductive (see [1], §6, no. 4). Hence go(0) is reductive as well. Set Go,/i = {g & 
Go \ 9 ■ h = h}. Then Go,h is the subgroup of G corresponding to go(0). Since a is a Z- 
graded reductive Lie algebra, gi(2) has a finite number of Go./i-orbits (see [29], §2.6), and 
therefore a unique dense one. Denote this orbit by U. Let u G gi(2); then the Go,h-orbit 
of u is dense if and only if it tangent space [go(0),M] is equal to gi(2). So we get that 
U = {u e gi(2) I [go(0),M] = gi(2)}. Now let e G gi(2) lie in a normal s[2-triple {h,e,f). 
Let K denote the Lie algebra spanned by h, e, /, i.e., K is isomorphic to sl2- Then o is a 
i^-module. From s[2-representation theory it follows that ade : go(0) gi(2) is surjective. 
Consider the orbit Go,/i ■ e C gi(2); its tangent space is [go(0), e] = gi(2). Therefore, e G f/. 
But then, since f/ is a single Go,/i-orbit, all m G f/ lie in a normal s[2-triple, {h,u, fu). So, 
if a given element in U does not lie in such a triple, then it follows that there is no normal 
5[2-triple containing h. □ 
This proposition immediately implies that the following algorithm is correct. 

Algorithm 8 Input: an admissible h G go. 
Output: True if h is normal, FALSE otherwise. 

1. Compute the spaces gi(2) and g_i(— 2). 

2. If h ^ [gi(2), g_i(-2)] then return False. 

3. By trying a few random elements find an e G gi(2) with [go(0), e] = gi(2). 

4. By solving a system of linear equations we decide whether there is an / G g_i(— 2) 
such that [h, e, /) is a normal sl2-triple. If such an / exists, then return True, 
otherwise return FALSE. 
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Remark. In the third step we need to find random elements. This can be done as follows. 
Let f2 be a finite set of integers (containing, say, all integers from to n, for some n > 0). 
Let Ui, . . . , Us be a basis of 0i(2). Then we choose randomly, independently, and uniformly 
s elements at E Q and form the element e = Yli=i (^i'^i- If ^ is large enough then e will lie 
the U from Proposition [7] with high probability; so we expect to find such an e within a 
few steps. 

The problem is how n should be chosen. In practice, usually a small n suffices. However, 
one can also proceed as follows. First we choose a small n. Then every time the random e 
does not lie in U we increase n. 

Also we remark that the correctness of the output does not depend on the way Vt is 
chosen. The output is always correct; only the running time depends on the choice of Vt. 

Remark. It is possible to dispense with the second step. However, practical experience 
has shown that the algorithm on the average becomes more efficient when we include it. 



3.3 Finding the nilpotent orbits 

Now we assume that the automorphism 9 is inner. Hence the Cartan subalgebra f)o of So 
is also a Cartan subalgebra of q. Since 0o is reductive we can write 0o = t © i^, where r is 
the centre of go and t = [0O)0o] is semisimple. Note that Gq = L x R, with L and R the 
subgroups corresponding to I and r respectively. The action of R on go is trivial; so the 
Go-action on go has the same orbits as the L-action. Set t)o = P)o H [; then [)q is a Cartan 
subalgebra of [. (Indeed, it is a maximal toral subalgebra of [, hence a Cartan subalgebra, 
as I is semisimple.) Let Wi = NL{i)Q)/CL{i)Q) be the Weyl group of [. This group acts on 
f)Q, and a closed Weyl chamber is a fundamental domain for this action. Let \E' be the root 
system of [ relative to i)^. Let H = {Pi, . . . , /3m} be a set of simple roots in Then the set 
of all h E i)Q with (3i{h) > for 1 < z < m is a fundamental Weyl chamber C[. The next 
lemma implies that we can assume that the candidates for normal /i € go lie in C[ © r. 

Lemma 9 Let (/ii,ei,/i) be a normal si2-triple. Then there is a unique normal 3i2-triple 
{h,e,f) that is Go-conjugate to {hi,ei, fi), and such that h = h' + u, with h' G C[ and 
u Ex. 

Proof. Observe that hi lies in a Cartan subalgebra of go, which is Go-conjugate to [)o- So 
by acting with an element of Go we find a conjugate normal s[2-triple {h2, 62, /2) such that 
^2 £ f)o- Write /i2 = ^3 + u, with G [)[, and m G r. Now h^ has a unique W^pconjugate 
/i4 lying in G[. So we find a normal s[2-triple (/i4 + u, 64, /4), conjugate to the original one. 
Moreover, a normal s[2-triple {h^ + v, 65, /s) with G G[, and f G r and 7^ /i4 cannot 
be Go-conjugate to (/ii,ei,/i), because two elements of t)o are L-conjugate if and only if 
they are VTrconjugate. □ 
Let $ be the root system of g with respect to ()o- Let W be the corresponding Weyl 
group. Then W acts on 1)q. It is well-known that the action of W on [)o can be realised 
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as follows. Let A = {ai, . . . , ai} be a basis of simple roots of $. Let h^^ be the unique 
element of ()o with aji^haj = (ctj, d^)- Then for the simple reflections Sq. we have Sa^h) = 
h — ai{h)hai- Let C C [)o be the fundamental Weyl chamber, i.e., C consists of all /i G [)o 
with ai{h) > for 1 < z < /. Since f)o is also a Cartan subalgebra of go, we have that \l/ 
is a root subsystem of $. This implies that Wi is a Weyl subgroup of W. We note that 
C( © r is the fundamantal Weyl chamber corresponding to the action of Wi on [)o. We 
choose the positive roots in \1/ such that they are also positive in $. Then C C C[ © r. 
Since Wi is a Weyl subgroup of W, by using the algorithm of Section \2A] we can compute 
a set {wi, . . . , Wm} of shortest length right coset representatives of Wi in W. 
Now we have the following algorithm. 

Algorithm 10 Input: an admissible h & C . 

Output: the list of all normal 5[2-triples (/i', e', /') with h' E Ci® r, and such that h' is 
VT-conjugate to h. 

1. Set Ti. = {wi{h) \ 1 < i < m}. (Here the Wi are the shortest length right coset 
representatives, as above.) 

2. Set T = 0. For each h E Ti we test whether h is normal using Algorithm [HI If this is 
the case then we find a normal s[2-triple containing h. We add it to T. 

3. Return T. 

Proposition 11 AlaorithmUU is correct. 

Proof. Note that by Lemma [21 the right coset representatives Wj are exactly the elements 
of W that map C into Ci © r. So TC is exactly the set of h' E CiQ) t that are ly-conjugate 
to h. From this set we take the normal elements. □ 
Remark. Note that it is straightforward to find an s[2-triple containing a normal element 
of f)o- In fact. Algorithm [8] finds one. 

The nilpotent G-orbits in g have been classified in terms of so-called weighted Dynkin 
diagrams (see [5], [6]). From such a diagram it is straightforward to find ah E C lying in an 
s[2-triple {h, e, /) (cf. [S]). This also implies that e is a representative of the corresponding 
nilpotent G-orbit. Let hi, . . . ,ht denote the elements of C so obtained. 

Lemma 12 Let h E i)o be normal. Then there is a unique hi such that h is W -conjugate 
to hi. 

Proof. The hi lie in sl2-triples {hi, Cj, fi). Let {h, e, /) be a normal s[2-triple. Then there is 
Ci such that e is G-conjugate to Cj. But that implies that the triples {h, e, /) and {hi, Ci, fi) 
are G-conjugate. Hence h and hi are G-conjugate. Now two elements of [)o are G-conjugate 
if and only if they are ly-conjugate. □ 
Now we apply Algorithm [10] to each of the hi, obtaining sets of Tj. Let T be their union. 
We claim that for each nilpotent Go-orbit in Qi the set T contains exactly one s[2-triple 
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{h, e, f) such that e is a representative of the orbit. Indeed, let e' G Qi be nilpotent, lying in 
a normal s[2-triple {h', e', /'). Then h' is Go-conjugate to an h" G f)o (as h' lies in a Cartan 
subalgebra of Qq, which is Go-conjugate to f)o). By the previous lemma h" is W^-conjugate 
to exactly one hi. But by Lemma [9l h" is Wo-conjugate to a unique /i G G[ © r. Hence 
Proposition [TT] implies that Tj contains an 5[2-triple (/i, e, /). We conclude that Tj contains 
at least one 512-triple {h, e, /) such that h is Go-conjugate to h'. If there were another one, 
{hi, ei, /i), then h and /ii would be Go-conjugate, and since they both lie in G[ © r, they 
would be equal by Lemma [91 

Remark. The algorithm also works when 6 is an outer automorphism. However, in that 
case Wi is no longer a subgroup of W; therefore in Algorithm [10], the set TC needs to be 
the intersection W ■ h (1 Ci (B x. The need to compute the entire orbit W ■ h makes the 
algorithm much less efficient. 

Remark. In [22j also a method based on computing the set Ti, as in Algorithm [TD] is 
developed. However, in order to decide whether a given h' G G( © r is normal, a more 
complicated procedure is used. 



4 Listing nilpotent orbits II 

In this section we describe an algorithm for finding representatives of nilpotent orbits using 
Vinberg's theory of carrier algebras. Again we divide the section into various subsections. 
We recall that g = ©^o^0« grading of g relative to the automorphism 9. Throughout 

we let f)o be a fixed Cartan subalgebra of go- By Wq we denote the Weyl group of the root 
system of go (relative to f)o)- 

4.1 Z-graded subalgebras 

We consider semisimple Z-graded subalgebras s of g, where Sk C gfemodm- The following 
terminology is used: 

• s is called regular if it is normalised by a Cartan subalgebra of go. 

• A regular s is called standard if it is normalised by P)o- 

• A regular s is called complete if it is not a Z-graded subalgebra of a Z-graded regular 
semisimple subalgebra of the same rank. 

• 5 is called locally flat if dimso = dimsi. 

Let s be a standard Z-graded semisimple subalgebra. Then s is contained in a unique 
complete standard Z-graded semisimple subalgebra, of the same rank ([30], Proposition 3). 
We call it the completion of s. 
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Let e' G 0i be nilpotent; then after raplacing e' by a Go-conjugate e we get that e hes 
in a complete standard locally flat Z-graded subalgebra s, called the carrier algebra of e. 
Moreover, e G 5i is in general position, which means that [sq, e] = 5i ([30], Theorem 2). 

Let e, e' G flo be nilpotent elements lying in carrier algebras 5, s'. As 5, s' are normalised 
by f)o5 their respective root systems are sets of weights of P)o; therefore, PFq acts on 

their elements. We have that e, e' are Go-conjugate if and only if are PVo-conjugate 

([30], Corollary to Theorem 2). 

Example 13 We consider the set up of Example [H Let e = 62,3 + 63^1. Let Sj for i = 
-2, -1, 0, 1, 2 respectivey be the spaces spanned by {e^a}, {63,2, 61,3}, {61,1-63,3, 62,2-63,3}, 
{^2,3, 63,1}, and {62,1}. Then 5 = s_2 © ■ ■ ■ © S2 is the carrier algebra of e. 

4.2 The completion 

Let 5 C be a standard semisimple Z-graded subalgebra. Set t = f)o H So; then t is a 
maximal torus in So- Indeed, if t could be enlarged to a bigger torus, then so could f)o, 
but the latter is a maximal torus in go. It follows that t is a Cartan subalgebra of So, and 
therefore of s (as s is Z-graded). Let 11 be a basis of simple roots of the root system of s 
with respect to t. Every root has a degree: deg(a) = if s^ C Sfc. We can choose 11 so 
that deg(a) > for all a G 11. Furthermore, if s is locally fiat then deg(Q;) G {0, 1} for 
a G n ([50], §4.2). 

Let ho E the such that a{ho) = deg(a) for G 11. Then 

Sk = {x E 5 \ [ho,x] = kx}. 

Therefore, ho is called a defining element of s. It is uniquely defined by the choice of t and 

n. 

Set 

3 = {-u G f)o I «(«) = for all a G 11}. (1) 

Note that since s is regular, the elements of 11 are weights of l)o, in particular we can view 
n as a subset of f)Q. Therefore, the definition of 3 makes sense. Furthermore, f)o = t©3. 
Indeed, t fl 3 = 0, dim 3 + dimt = dim [)q. 
Now define the Z-graded subalgebra u by 

Uk = {x e Qk mod m I [d,x] = and [ho,x] = kx}. 

Proposition 14 u = s' © 3, where s' is a standard complete Z-graded subalgebra of the 
same rank as s, and containing s. In other words, s' is the completion of s. 

Proof, (cf. [30j, Proposition 2). Note that [3,s] = 0, and hence s C u. Furthermore, t©3 
is also a Cartan subalgebra of Uo. Hence the rank of u is equal to the rank of go- By [30] . 
Lemma 2, u is reductive. Hence 3 is the centre of u. So s' = [u, u] is semisimple, contains 
s, and is of the same rank as s. Moreover, by [3Qj, Proposition 1, s' is complete. Finally, 
as P)o C u, also s' is standard. □ 
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Now for the remainder of this section we suppose that the automorphism 6 is inner. 
Then ()o is also a Cartan subalgebra of g. Therefore, s is spanned by root spaces of g. By 
$ we denote the root system of g with respect to [)o. Furthermore, for a G $, go, will be 
the corresponding root space. 

Set 

^0 = {« e $ I ga C 00 and a{ho) = 0, a{^) = 0}, 

^'i = {a e $ I ga C 01 and a{ho) = 1, 0(3) = 0}. (2) 

Let s' be the completion of s. Then from Proposition [T3] it follows that Sq is the sum of t 
and the 0^ for a G \E'o- Also, s'^^ is the sum of the go, for a G \E'i. So s' is locally flat if and 
only if |n| + |\E'o| = |^i|- Moreover, Hq is also a defining element of s'. 

We summarise the above in the following algorithm. The input will be a basis of simple 
roots n C $ of the root system of a standard semisimple Z-graded subalgebra s. We assume 
that deg(a) G {0, 1} for a G 11 and that we have given a decomposition 11 = Ho U IIi, 
where Hi contains the roots of degree i. The output will be a defining element of the 
completion s' if s' is locally flat. Otherwise the output will be fail. The algorithm takes 
the following steps: 

1. Compute a basis of i = f)o HSq. (This can for instance be done by taking root vectors 
Xa G ga, X-a G g_Q, for a G 11; then t will be spanned by the [xq,,x_q,].) 

2. Compute Hq G t, with a(/io) = deg(a) for a G 11. 

3. Compute a basis of 3 C f)o, where 3 is as in ([T]). 

4. Compute \i/o and ^'i, as in ([2]). 

5. If |n| + |\l'o| 7^ l^il then return fail. Otherwise return Hq. 

Example 15 Let g, s be as in Examples [H [131 After some small calculations we get 
ho = — ei^i + ei_2 and 3 is spanned by Ci^i + 62,2 + 63,3 — 364^4. Using this it is straightforward 
to see that the completion of s is s itself. 

4.3 Obtaining the nilpotent orbits 

On the basis of what was said in the previous subsections we describe an algorithm for 
listing the nilpotent Go-oi'bits in gi. Also in this subsection we assume the automorphism 
9 to be inner. Let $0 (respectively, $1) be the set of roots a of g such that ga C go 
(respectively, C 0i). We note that $0 is a root sub-system of $. Let Aq be a basis of 
simple roots of $o- 

In the first step in our algorithm we obtain a set V of vr-systems (cf. Section 12. 3p 
contained in $0 U $1. This set is required to contain, up to Wo-conjugacy, bases of the 
root systems of all locally fiat, standard, semisimple Z-graded subalgebras of 0. However, 
the set may be bigger than necessary, i.e., it may include M/Q-conjugate vr-systems, and it 
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may have vr-systems that are bases of root systems of Z-graded subalgebras that are not 
locally fiat. 

First we note the following. Let 11 = Ho U Hi be a basis of the root system of a standard 
semisimple Z-graded subalgebra of g. Here Ilj C contains the roots of degree i. Then 
Ho is a TT-system in $o- This means that the following steps will produce a suitable list V: 

1. Use the algorithm in Section 12.31 to get the set Vq of vr-systems in $0) up to Wq- 
conjugacy. 

2. For all Ho G Vq we find all maximal Hi C $i such that Ho U Hi is a vr-system. This 
yields a set V'l of vr-systems in <l>o U $i. 

3. From V'l we erase PFo-conjugate copies (see Section [2l2l) . to get a set V". 

4. Finally, set "P = 0. For all Uq U Hi in we add all Uq U U[ to V, where U[ runs 
through the subsets of IIi. 

In the second step of our algorithm we find a list of normal s[2-triples corresponding 
to the nilpotent Go-orbits in gi. This step consists in running through the list V, where 
for every IIo U IIi we do the following. First we apply the algorithm of Section 14.21 If 
the result is fail, then we go to the next element of V. Otherwise the result is a defining 
element ho G f)o of the completion s' of the Z-graded semisimple subalgebra of g with root 
system spanned by Hq U Hi. Furthermore, s' is locally fiat. As outlined in Section 14111 s' 
corresponds to a nilpotent Go-orbit in Qi. We obtain this orbit as follows. We observe that 
h = 2hQ is such that there is a normal 5l2-triple {h, e, /), where e is a representative of the 
orbit we are after ([30], §4.2). Using the notation of Section [3731 we compute the unique 
h E CiQ) t that is PVo-conjugate to h. There is a normal s[2-triple {h, e, /), where e lies in 
the orbit we are interested in. If the element h has not occurred before, then we find such 
an s[2-triple, and add it to our list. Otherwise, we go to the next element of V. 
Remark. On the algorithms of this section some variations are possible. Firstly, in 
the set V we can erase all PVo-conjugate copies of vr-systems, not just of the maximal 
ones. This yields a much smaller set, which makes the second step easier to execute. 
However, practical experiences show that this leads to a less efficient algorithm (i.e., there 
is much more work involved in checking VFo-conjugacy for all pairs of elements of V, than 
in computing the elements h G G[ © c). The advantage of only checking the maximal 
vr-systems for Wo-conjugacy, is that with each instance also a lot of subsystems are shown 
to be Wo-conjugate. 

Secondly, the algorithm of Section 14.21 could just return the ho, and not worry about 
flatness. In that case, if the corresponding graded subalgebra is not locally flat, then 2hQ 
may not correspond to a nilpotent orbit. This can then be checked using the methods of 
Section [31 However, practical experience shows that this leads to a less efficient algorithm 
as well. 
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Remark. It is straightforward to extend the algorithm to outer autmorphisms. Indeed, 
in that case, $o is the same as before, wheras $i becomes the set of weights of f)o in Qi. 



5 Practical experiences 

The algorithms of this paper have been implemented in the language of GAP4. In this 
section we report on their running times on some sample inputs. As input we have taken 
a few N-regular automorphisms of the Lie algebra of type Ej (see Section [6]). The running 
times are displayed in Table [H The set H is the union of all sets, denoted with the same 
letter, that occur in the first step of Algorithm [10] (recall that this algorithm is called 
several times). This means that Ti. is the set of all admissible elements of f)o that are tested 
for normality. Secondly, the set V is as in Section 14.31 





Method I 


Method II 






e 






Wi\W\ 


\n\ 


time 


\V\ 


time 


dim 00 


dim 01 


# orbits 






72 


721 


41 




oo 


63 


70 


94 






672 


4627 


83 


4227 


529 


43 


45 


75 


4 


4032 


22939 


475 


4014 


165 


33 


35 


113 






10080 


52109 


1650 


2494 


31 


27 


27 


82 


6 


40320 




oo 


4302 


50 


21 


24 


233 



Table 1: Running times of the algorithms, with input an N-regular automorphism 6 of the 
Lie algebra of type Ef. The first column has the order of 6. The next three columns display 
data relative to the algorithm of Section [3l the second column has the index of W[ in W, 
the third column has the size of Ti, and the fourth column lists the time taken. The next 
two columns display data relative to the algorithm of Section HJ the fifth column has the 
size of the set V, and the sixth column the total time taken. The last three columns list, 
respectively, the dimensions of 0o and 0i and the total number of nilpotent orbits found. 
All running times are in seconds. An oo indicates that the computation did not terminate 
within one hour. 

From the Table we see that Method I behaves well when the dimension of the space 0o 
is big, because then the index of W[ in W is small, which results in far fewer admissible 
elements of t)o that need to be checked for normality. On the contrary. Method II behaves 
badly in that the set V gets too big. However, this method gets quickly better 

when the dimensions of the spaces 0o, 0i decrease. So indeed the two methods complement 
each other. 
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6 N-regular automorphisms 



Let m > 2 be an integer. By results of Antonyan and Panyushev (see [25j) there is a 
unique (up to conjugacy) inner automorphism of order m such that Qi contains a regular 
nilpotent element. For outer automorphisms a similar statement holds (see [25]). For the 
Lie algebra of type Eq this means that if there are outer automorphisms of order m, then 
there is exactly one (up to conjugacy) such that gi contains a regular nilpotent element. 

Automorphisms such that Qi contains a regular nilpotent elements are called N-regular. 
These N-regular automorphisms have a number of special properties (see [25]), but it is 
not immediately obvious which inner automorphisms are N-regular. 

By using the algorithms of the previous sections we can find the N-regular automor- 
phisms of order m. Indeed, we can list representatives of all conjugacy classes of inner 
automorphisms of order m. For each element of the list we obtain the nilpotent Go-orbits 
in 01. In particular, we find representatives of these orbits. By computing their weighted 
Dynkin diagrams we check whether a regular nilpotent element occurs among them (we 
refer to [15] for an algorithm to compute weighted Dynkin diagrams). 

In Tables E], HJ [5], [6], and [7] we give the Kac diagrams of the N-regular inner automor- 
phisms of the Lie algebras of exceptional type, of orders between 2 and h — 1 (where h is 
the Coxeter number). Furthermore, Table [3] contains the Kac-diagrams of N-regular outer 
automorphisms of the Lie algebra of type Eq. These Kac diagrams are extended Dynkin 
diagrams with labels that define the automorphism. It turns out that in our tables only 
the labels and 1 occur. Therefore, we give the Kac diagram by colouring the nodes of 
the extended Dynkin diagram. A black node has label 1, a non-black node has label 0. 

The contents of the tables are as follows: the first column contains the order of the 
N-regular automorphism, and the second column its Kac diagram. The third column lists 
the number of nilpotent orbits in Qi. Let Af denote the variety of all nilpotent elements in 
01. It is known ([19], [30]) that Af splits in irreducible components J\fi, where 1 < i < r. 
The component A/i is the closure of an orbit GoCi, where G 0i is nilpotent, such that 
the orbit GoCi is of maximal possible dimension. Therefore, all components are of the 
same dimension, and we can compute the Cj as follows. First we list representatives of all 
nilpotent orbits in gi. Let e be such a representative; then the dimension of its orbit is equal 
to the dimension of [go, e]. Hence, by simple linear algebra, we can compute the dimension 
of each orbit. The Cj are then the representatives of the orbits of maximal dimension. In 
particular, we get the number of irreducible components, and their dimension. That is the 
content of the fourth and fifth column. The last column displays the rank of gi, that is 
the codimension of a nilpotent orbit of maximal dimension. (It is also the dimension of a 
Cartan subspace, cf. [29]). 

When the order of 6 is 2, then all nilpotent orbits in gi of maximal dimension are 
conjugated in the group G ([19], Theorem 6). A similar statement fails to hold for larger 
orders. Following [21] , we say that an A^-regular automorphism for which all nilpotent Go- 
orbits in gi of maximal dimension lie in the same G-orbit, is very N-regular. We checked 
which 6 are very N-regular by computing the weighted Dynkin diagrams of the orbits Gci, 
where GoCi is a nilpotent Go-orbit in gi of maximal dimension. If 6 is not very N-regular 
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then the number in the fourth column has an added *. The situation where M has more 
than one component, and 9 is very N-regular, appears to occur very rarely. 
Remark. There is a also a more direct method for constructing the N-regular inner 
automorphisms of order m (cf. [Ij). Let e G g be a regular nilpotent element lying in an 
5[2-triple (/i, e, /). Since e is even, all eigenvalues of adg/i are even integers. Let 



00(2fc) 







be the grading of g into adg/i-eiganspaces. For < z < m — 1 let be the sum of all 



g(2/c) such that k mod m = i. Then g 



^m—l 



Qi is the grading of g corresponding to an 



A^-regular automorphism of order m. It is inner as go contains a Cartan subalgebra of g. 
Some work is still needed to obtain the Kac diagram from the grading. 

Table 2: N-regular inner automorphisms of Eq. 
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N-regular automorphisms of Eq. 
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Table 3: N-regular outer automorphisms of £^6- 
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Table 4: N-regular inner automorphisms of E^. 
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N-regular automorphisms of E^. 
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Table 5: N- regular inner automorphisms of Eg. 
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N-regular automorphisms of Eg. 
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N-regular automorphisms of Eg. 
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Table 6: N-regular inner automorphisms of F4. 



order 



Kac diagram 



# orbits 



# components 



dim 



rank 



2 
3 
4 

5 
6 
7 
8 
9 

10 
11 



o 



o 



o 



o 



-o 



-o 



-o 



-o 



-o 



-0= 



^^^^ 



3> 



26 
19 
29 
15 

35 
13 
30 
19 
19 
15 



1 
1 

3* 
1 

6* 
1 

4* 
1 
1 
1 



24 

16 

12 

11 

8 

7 

6 

6 

6 

5 



4 
2 
2 

2 

1 






Table 7: N-regular inner automorphisms of G2- 



order 



2 
3 

4 
5 



Kac diagram 



O 



-CE 



O 



# orbits 



5 
6 

4 
3 



# components 



1 

2* 
1 
1 



dim 



6 
4 

4 
3 



rank 



2 
1 
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